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Abstract 

It is shown that the CKMT model for the nucleon structure function F 2 , taken 
as the initial condition for the NLO evolution equations in perturbative QCD, 
provides a good description of the HERA data when presented in the form 
of the logarithmic slopes of F 2 vs x and Q 2 (Caldwell-plot), in the whole 
available kinematic ranges. Also the results obtained for the behavior of the 
gluon component of a nucleon are presented. 
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1 The CKMT model 

The CKMT model |1| for the parametrization of the nucleon structure function F 2 
is a theoretical model based on Regge theory which provides a consistent formulation of 
this function in the region of low Q 2 , and describes the experimental data on F2 in that 
region. 

The CKMT model |T| proposes for the nucleon structure functions 

F 2 (x, Q 2 ) = F s (x, Q 2 ) + F NS (x, Q 2 ), (1) 

the following parametrization of its two terms in the region of small and moderate Q 2 . 
For the singlet term, corresponding to the Pomeron contribution: 

/ ra \ i+a(q 2 ) 

F s (x, Q 2 )=A- x-^ ■ (1 - x)<^ ■ [^-\ , (2) 

where the x— >0 behavior is determined by an effective intercept of the Pomeron, A, which 
takes into account Pomeron cuts and, therefore (and this is one of the main points of the 
model), it depends on Q 2 . This dependence was parametrized as : 

Thus, for low values of Q 2 (large cuts), A is close to the effective value found from analysis 
of hadronic total cross-sections (A~0.08), while for high values of Q 2 (small cuts), A takes 
the bare Pomeron value, A~0.2-0.25. The parametrization for the non-singlet term, which 
corresponds to the secondary reggeon (f, A 2 ) contribution, is: 



F NS (x, Q 2 )=B- x 1 ^ ■ (1 - x)"^ 2 ) • ( ) , (4) 



Q 2 + b / 



where the x— >0 behavior is determined by the secondary reggeon intercept atR, which 
is in the range a^=0.4-0.5. The valence quark contribution can be separated into the 
contribution of the u (B u ) and d (Bj) valence quarks, the normalization condition for 
valence quarks fixes both contributions at one given value of Q 2 (we use Q 2 = 2.GeV 2 
in our calculations). For both the singlet and the non-singlet terms, the behavior when 
x— >1 is controlled by n(Q 2 ), with n(Q 2 ) being 

so that, for Q 2 =0, the valence quark distributions have the same power, given by Regge 
intercepts, as in the Quark Gluon String Model or in the Dual Parton Model ||, 
n(0)=OR(0)— aAr(0)~ 3/2, and the behavior of n{Q 2 ) for large Q 2 is taken to coincide 
with dimensional counting rules. 

The total cross-section for real (Q 2 =0) photons can be obtained from the structure 
function F2 using the following relation: 



q 2 •F 2 {x,Q ) 



(6) 
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The proper F 2 (x, Q 2 ) r ^Q 2 behavior when Q 2 — >0, is fulfilled in the model due to the last 
factors in equations |2| and (|. Thus, the a^iv) has the following form in the CKMT 
model: 

o™{v) = A7r 2 a EM ■ (A ■ a- l - Ao • (2mi/) Ao + (B u + B d ) ■ b' aR ■ {2mv) aR - 1 ) . (7) 

The parameters were determined flp from a joint fit of the <t*°* data and the NMC 
data [U on the proton structure function in the region lGeV 2 < Q 2 < 5GeV 2 , and a very 
good description of the experimantal data available was obtained. 

The next step in this approach is to introduce the QCD evolution in the partonic 
distributions of the CKMT model and thus to determine the structure functions at higher 
values of Q 2 . For this, the evolution equation in two loops in the MS scheme with 
A = 200. MeV was used 0. 

The results obtained by taking into account the QCD evolution in this way are [0 in 
a very good agreement with the experimental data on F 2 (x, Q 2 ) at high values of Q 2 . 

When the publication of the data || || on F 2 from HERA at low and moderate Q 2 
provided the opportunity to include in the fit of the parameters of the model experimental 
points from the kinematical region where the CKMT parametrization should give a good 
description without need of any perturbative QCD evolution, one proceeded JF| to add 
these new data on F 2 from HI and ZEUS at low and moderate Q 2 , to those from NMC 
H and E665 |1 collaborations, and to data M on cross-sections for real photoproduction, 
into a global fit which allowed the test of the model in wider regions of x and Q 2 . One 
took as initial condition for the values of the different parameters those obtained in the 
previous fit [0], and although the quality of the fit is not very sensitive to small changes in 
the values of the parameters, the best fit has been found for the values of the parameters 
given in Table(l|. 

Table 1: Values of the parameters in the CKMT model obtained in the fit of F 2 when 
also the low Q 2 HERA data are included. All dimensional parameters are given in GeV 2 . 
The valence counting rules provide the following values of B u and B^, for the proton case, 
when fixing their normalization at Ql=2.GeV 2 : i? u =1.1555, £^=0.1722. 



CKMT model 


values of the parameters 


A 


0.1301 


a 


0.2628 


A 


0.09663 


Ax 


1.9533 


A 2 


1.1606 


c 


3.5489 (fixed) 


b 


0.3840 


an. 


0.4150 (fixed) 



The quality of the description provided by the CKMT model of all the experimental 
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data on cr*° and F 2 , and, in particular, of the the new experimental data from HERA is 



very high, with a value of \ 2 /d.o.f. for the global fit, x 2 /d-o.f. =106. 95/167, where the 
statistical and systematic errors have been treated in quadrature, and where the relative 
normalization among all the experimental data sets has been taken equal to 1. 

Thus, by taking into account the general features of the CKMT model described 
above, we use the CKMT model to describe the experimental data in the region of low Q 2 
(0 < Q 2 < Qq), and then we take this parametrization as the initial condition at to be 
used in the NLO QCD evolution equation to obtain F 2 at values of Q 2 higher than Qq. 
In order to determine the distributions of gluons in a nucleon the CKMT model assumes 
[ ij that the only difference between distributions of sea-quarks and gluons is in the x — > 1 



behavior. Following |T(J we write it in the form 

xg{x,Q 2 ) =Gxq{x,Q 2 )/{l-x), (8) 

where xq(x, Q 2 ) is proportional to the expression in equation ||. The constant G is deter- 
mined from the energy-momentum conservation sum rule. 

We have performed our calculations for two different values of = 2.GeV 2 and 
Ql = A.GeV 2 . We also show our results in the shape of both the dF 2 /dlnQ 2 and the 
dlnF2 j dlnil j x) slopes in order to compare with the experimental data when given in the 
so-called Caldwell-plot. This approach provides a smooth transition from the region of 
small Q 2 , which is governed by the physics of Regge theory, to a region of large Q 2 , where 
the effects of QCD-evolution are important. 

The way we proceed to calculate F 2 , and its logarithmic derivatives dF 2 / 'dlnQ 2 , and 
dlnF2/dln(l/x), is the following (see Appendix A for all the technical details on how the 
NLO QCD evolution has been performed): 

• In the region < Q 2 < Qq we use the pure CKMT model for F 2 . 



• For Qq < Q < charm threshold |TT|, we make the QCD evolution of F 2 at NLO 
in the MS scheme for a number of flavours rif = 3, and we take as the starting 
parametrization that given by the CKMT model. For Ql we have used in this 
calculation two different values: Ql = 2.GeV 2 , and Ql = A.GeV 2 . 

• When charm threshold < Q 2 < Q 2 = 50. GeV 2 , also the QCD evolution of F 2 is 
implemented at NLO in the MS scheme for a number of flavours rif = 3, using 
the parton distribution functions for the u, d, s quarks, and by including the charm 
contribution via photon-gluon fusion. 

• For values of Q 2 > Q 2 , QCD evolution is computed at NLO in the MS scheme, 
but now with a number of flavours rif = 4, and by using the parton distribution 
functions for the u, d, s, and c quarks. 

One has to note that in the treatment of the charm contribution we have followed 
reference 1T2 . 



2 Results 

The results we have obtained are presented in figures 1 to 9. 
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Figure 1 shows F 2 (x, Q 2 ) as a function of x for several values of Q 2 , from Q 2 = 0.6GeV 2 
to Q 2 = 17.GeV 2 . The dotted lines correspond to the pure CKMT model without any 
perturbative evolution, while the full lines run for the evoluted CKMT parametrization. 
When for a given value of Q 2 two full lines are depicted, the bold (solid) one has been 
obtained by taking the starting point for the QCD evolution as Ql = 2.GeV 2 (Qq = 
A.GeV 2 ). Experimental points in this figure are from E665 HI flT3" |, and ZEUS [|TJj] 
collaborations. 

In Figures 2. a and 2.b, we present the comparison of the pure CKMT parametrization 
of F 2 with the low Q 2 data of E665, ZEUS-BPC95, and ZEUS-BPT97, as compiled in 



and fl6fl . One sees that the agreement between the CKMT model and the experimental 
data in this region of low Q 2 is good. 

In Figure 3 (Caldwell-plot), the slope dF2/dlnQ 2 is shown as a function of x, and 
compared with the a + blnQ 2 fit to the ZEUS F 2 data in bins of x. This plot was 
considered as an evidence for a transition from hard to soft regime of QCD in the region 
of Q 2 ~ 5.GeV 2 (see for example |T7|]). This question has been studied theoretically in 
references [[18|, |19[ . Figure 3 shows that the CKMT model is in a good agreement with 
experimental points in the whole region of x and Q 2 . One problem with the presentation 
of the data in Figure 3 is a strong correlation between x and Q 2 values for the data points. 
It follows from the formulas of CKMT model for dF 2 (x, Q 2 )/dlnQ 2 given in Appendix B 
that for a fixed value of Q 2 this quantity monotonically increases as x — > 0. The existence 
of a maximun of dF 2 (x,Q 2 )/dlnQ 2 in Figure 3 is related to the correlation between Q 2 
and x in the region of small x (or Q 2 ). The same conclusion was achieved in reference 



Tq| , and recently confirmed by experimental data JT 

Figures 4 and 5 show the slope dlnFlj dln{\j x) as a function of Q 2 compared to the 
fits F 2 = Ax~ Ae ff of the the ZEUS |L4| and HI [[yj data, respectively. In Figure 4, as 
the x range of the BPC95 data is restricted, also the E665 H] data were included in WA 



and are now also taken into account. This slope is sometimes interpreted as the A e ff of 
the Pomeron exchange, A e // = dlnFlj dln{\j x). Let us note that in our approach A e ff 
for Q 2 > Qq can not be interpreted as an effective Pomeron intercept, because the QCD 
evolution leads to a substantial increase of A e jj as Q 2 increases. On the other hand this 
effect should decrease as x — ► 0. 

In the experimental fits, each Q 2 bin corresponds to a average value of x, < x >, 
calculated from the mean value of ln(l/x) weighted by the statistical errors of the cor- 
responding F 2 values in that bin. Even though we can proceed as in the experimental 
fits, and we get a very good agreement with the data, since the estimation of < x > is 
in some sense artificial and arbitrary, and it introduces unphysical wiggles when drawing 
one full line connecting the different bins, we made for all the Q 2 bins in this figures the 
choice of the smallest x in the data, instead of considering a different < x > for each Q 2 . 
This choice is based on the fact that the ansatz Ag// = dlnF2 / dln(l j x) is actually valid 
for small x, and it results in a smooth curve except for the jump in the region around 
Q 2 ~ 50.GeV 2 , where the evolution procedure changes (again, see Appendix A for more 
details). 

Since the structure function F 2 in the region of low x is determined at large extent by 
the gluon component, we present our prediction for the behavior of this gluon component. 
Thus, Figure 6 shows the gluon density distribution as a function of Q 2 calculated by 
performing the NLO QCD evolution of the CKMT model, and its comparison with the 
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HI Collaboration data in reference and Figure 7 represents the gluon densities at 
/i 2 = 25.GeV 2 as a function of x calculated by evoluting the CKMT model at NLO in the 
QCD evolution, and compared to those determined from HI DIS and photoproduction 
data. Experimental data on D* meson cross-section measurements are from references 
16. p0[| . Figure 8 shows the behavior of xg(x,[i 2 ) at fi 2 = 200 GeV 2 as a function of Q 2 



to be compared with the Hl-dijets results [IB, [HJ. Finally, Figure 9 shows the prediction 
of the CKMT model for xg(x, Q 2 ) as a function of x at the values of Q 2 measured both 
by HI and ZEUS collaborations. 

A satisfactory agreement with the experiment is obtained in the whole ranges of x 
and Q 2 where experimental data are available, showing that the experimental behavior 
of F 2 , its logarithmic slopes, and its gluon component can be described by using as initial 
condition for the QCD evolution equation a model of F 2 where the shadowing effects 
which are important at low values of Q 2 are included, like the CKMT model. 



3 Conclusions 

The CKMT model for the parametrization of the nucleon structure functions provides 
a very good description of all the available experimental data on F 2 (x,Q 2 ) at low and 
moderate Q 2 , including the recent small- £ HERA points. 

An important ingredient of the model is the dependence of an effective intercept of 
the Pomeron on Q 2 . It has been shown recently |22] that such a behavior is naturally 
reproduced in a broad class of models based on reggeon calculus, which describes simul- 
taneously the structure function F 2 and the diffractive production by virtual photons. 

Use of the CKMT model as the initial condition for the QCD-evolution equations in 
the region of Q 2 = 2. 4- 5. GeV 2 leads to a good description of all available data in a 
broad region of Q 2 , including the logarithmic slopes of the structure function F 2 (x, Q 2 ), 
dF 2 (x, Q 2 )/dlnQ 2 and dlnF 2 (x, Q 2 ) / dln(l / x) . Thus an unified description of the data on 
F 2 for all values of Q 2 is achieved. 
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Appendix A - NLO QCD evolution of F 2 (x,Q 2 ) 



For the reader convenience we present here some technical remarks concerning the 
NLO QCD calculation of F 2 (x, Q 2 ). 

For sufficiently large Q 2 > 1 GeV 2 , the structure function F 2 (x, Q 2 ) can be expressed 
by perturbative parton distributions. In leading order (LO) perturbation theory, the 
expression is given as 



x 



F 2 (x, Q 2 )=xY: e\{q{x, Q 2 ) + q{x, Q 2 )}, 



(9) 



where q and q denote the quark and anti-quark distribution functions, e 2 the square of 
the quark electric charge, and the sum runs over all quark flavors included [|ll|]. On 
the other hand, with F 2 (x,Q 2 ) given in eqs. (|T]-[5|), and making reasonable assumptions 
concerning the flavor structure of the QCD-sea, one can extract from F 2 (x, Q 2 ) the differ- 
ent parton distribution functions, including that of the gluon component |T|]. Generally, 
the calculation of F 2 (x,Q 2 ) at Q 2 > 1 GeV 2 requires a Q 2 -evolution a la DGLAP f2~3fl . 
The procedure consists in the solution of the LO-DGLAP equations for the parton dis- 
tribution functions using reasonable initial distributions at a starting value Q 2 = 
(1 GeV 2 < Ql < 5 GeV 2 ). Using eq.(H), the resulting quark distributions at Q 2 can be 
recombined to F 2 at this virtuality. 

By the evolution of the CKMT-model we mean the application of this procedure to 
the model discussed in this paper. As mentioned above, the CKMT-model of F 2 (x, Q 2 ) is 
valid within < Q 2 < 5 GeV 2 . Due to the good aggreement with experimental data the 



parton distributions extracted from K 



CKMT 



at a Qq in the range given above seem to be 



reasonable initial distributions for an evolution to higher Q 2 . 

In next to leading order (NLO), the relation between F 2 (x, Q 2 ) and the parton distri- 
bution functions is more complicated and depends on the renormalization scheme. The 
calculations presented here are performed in the MS-scheme f24|. In this context, the 
structure function is given by [O as 



±F 2 (x, Q 2 ) = £ e 2 {?(*, Q 2 ) + q(x, Q 2 ) + ^p- [C q>2 * (q + q) + 2 • C g , 2 * g]) , (10) 



where q, q and g are the NLO quark, anti-quark and gluon distribution functions, respec- 
tively. a s denotes the strong coupling constant in NLO. The convolutions C * q and C * g 
are defined as 



C*q= f'^C (*) q(y,Q 2 ). 

J* y \yj 

The Wilson coefficients C qt g t2 (z) are given by 



CM = - 



1 + z 2 / 1 - z 3\ \ , 
In + 9 + 5 2 

1 - z \ z 4/4 



C,M = 2 



[ z 2 + (l- z ) 2 )]n 



1 - z 



- l + Sz(l - z) 



'121 



S 



Here, the integral over a [-^-distribution is denned as described in P5" 

-i dy 



C + * q 



-c 



x 

—< I- 



x y 
dy 



q(y,Q 2 ) 



x y 



x 



q(y,Q 2 )--q(x,Q 2 ) 

y 



q(x,Q 2 ) / dyC{y). (13) 



There are alternative renormalization schemes as, for instance, the DI S-scheme (II 
Here, the form of eq. (ffi) is kept for NLO also, i.e. 



-F 2 (x,Q 2 

x 



X 



J2 e q{<lDis(z, Q 2 ) + q D is(x, Q 2 )}. 



(14) 



The relation between the MS- and the D/S-distributions is given by 

• q ,2* V +C g ,2*9 +0{a 2 s ). 

Z.7I 

fn2\ r 

9dis(x,Q 2 ) = g(x,Q 2 ] 



{( 1dis (x, Q 2 ) 



V(,,g 2) + ^ 

a s (Q 2 ] 



(15) 



2tt 



J2 C ^*(q + q) + 2f-C g , 2 *g 



+ 0(ai 



The parameter / denotes the number of active flavors in the sea. 

Our procedure to extract the parton-distributions from p^ KMT i s based on the LO- 
formula eq. (0). Therefore, in NLO, we extract the -D/S-distributions. Now, the task 
is to calculate the MS-distributions at Q 2 = Q\ ■ This can be done using a first order 
approximation in a s (Q 2 )/27r: 



V (x,Q 



g{x,Ql 



q DIS [X, Qo) 

2n 



G q> 2* q DIS +G g ,2 * gDIS 



(16) 



9dis{x,QI 



C i? * (Qdis + Qdis) + 2/ • C 9i2 * gDis 



i q 



In summary, the Q 2 -evolution of F^ KMT works as follows: 

1. One chooses an appropriate value Q 2 = Qq > 1 GeV 2 as a starting point for the 
evolution. These are Q\ = 1 GeV 2 or Q\ = 4 GeV 2 in our calculations. 

2. At Q 2 = Qq, one extracts the NLO parton distribution functions from F^ KMT . 
The relation between these parton distributions and the structure function is given 
by eq.(Jl4]) which is formally the same as eq.(f|) in LO. So the resulting parton 
distributions are the D/S'-functions, i.e. qDis(x,Ql), q~Dis( x ,Qo) an d 9dis(x,Qq). 

3. Using eq.fllED one calculates the MS'-distributions q(x, Qq), q(x, Qq) and g(x, Qq). 



4. These MS'-functions serve as initial distributions in a numerical procedure to solve 
the NLO-DGLAP-equations in the MS'-scheme for a certain value Q 2 > Qq. The 
result are the evoluted MS-parton distributions q(x,Q 2 ), q(x,Q 2 ) and g(x,Q 2 ). 



5. Finally, using eq.flTOD, the structure function F 2 G (x, Q ) can be recalculated 
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The charm production is of particular interest. Following refs. [ |TT| , |T2| , the assumption 
of a "massless" charm quark produced above the threshold Q 2 = 4m 2 (m 2 - charm quark 
mass) via the usual DGLAP-evolution is not realistic. This procedure is useful in the 
range of high Q 2 3> Q 2 only. In the intermediate region Q 2 < Q 2 < Q 2 = 50 GeV 2 , the 
charm is treated via a photon-gluon fusion process. The corresponding contribution to 
the structure function is defined as 



\F;(x,Q\ml) = H^p-fJ^ ■C'J^)-g( V ,^), (17) 



where /i 2 = 4m 2 , a = 1 + 4m 2 /Q 2 , and, the coefficient C c g2 {Z,K) is given by 

C c gt2 (Z,R) = ^{[Z 2 + (l-Z) 2 + 4ZR(l-3Z)-8Z 2 R 2 ]\n^±^ 

+ V[-1 + 8Z(1-Z)-4ZR(1-Z)]}, (18) 

with V 2 = 1 — ARZ/ (1 — Z). So, F 2 in total is given by eq.([T0|) where the sum runs over 
q = u,d,s plus eq. (|TTD . The contributions of the botom and top quarks are neglected 
here. Precisely, the charm threshold is defined as discussed in refs. fill, IHfl > 



W 2 = Q 2 (l/x - l)>Q 2 =4m 2 . (19) 
The Q 2 -dependence of F 2 can be summarized as follows: 

i) Q 2 < Q 2 : In the low Q 2 region, F 2 is calculated as given in the pure CKMT-model, 

eqs.(|-0). 

ii) Ql < Q 2 < Q 2 : Below the charm threshold F 2 is calculated using eq. ([T0|) from NLO 

QCD evoluted parton distributions in the MS'-scheme. The number of flavors is 
/ = 3 (u,d,s). 

iii) Q 2 < Q 2 < Q = 50 GeV 2 : Above the charm threshold F 2 is determined from eqs. flTU|) 

and ([17]) . Note that the number of flavors active in the evolution is again / = 3 
(u,d,s). However, / = 4 after the charm is produced. This is important for the 
calculation of a s . 

iv) Q 2 > Q 2 = 50 GeV 2 : In the high Q 2 -region, F 2 is given by eq . (|10|) . The charm is 

produced as "massless" quark in the evolution process. Generally, the number of 
flavors is / = 4. 

The threshold Q 2 where the charm production in the evolution process is more important 
than the photon-gluon fusion is discussed in detail in |12| . The value of 50 GeV 2 is chosen 
to guarantee a transition as smooth as possible. This method works better for x — > 
than for x — > 1. This explains the small wiggles in some of the figures at Q 2 = 50 GeV 2 . 
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Appendix B — The slopes of i^^Q 2 ) 



For low Q 2 , the "pure" CKMT-model is used, i.e. the one defined in eqs.(|l|-|5]). Here, 
the calculation of the slopes as dF 2 (x, Q 2 )/d\nQ 2 and d\nF 2 (x, Q 2 )/d\n (1/x) = A e // is 
straightforward. Considering x and Q 2 as independent variables one gets 

dF t n § 2) = F s (x,Q 2 ) [^(A(Q 2 )-A )/n^ 



^ c ^m-l)ln(l-x) + a -^pl] (20) 



+ F NS (x,Q 2 ) (n(Q 2 ) - |) ln(l - x) + 



b Qfl(0) l 
=+6 



which in the limit Q 2 — > takes the form 

dF dl%$ 2) ~ (l + A )F s (x,Q 2 )+a R (0)F NS (x,Q 2 ). (21) 

Also, if one considers the case when W is fixed one can take x ~ C ■ Q 2 , and then, up 
to constant factors, one gets: 



-A(Q 2 ) + ^ (n{Q 2 ) - |) ln(l - Q 2 

Q 2 n(Q 2 ) i g (l+A(Q 2 )) i 
1-Q 2 ' Q 2 +a J 

+ i^ 5 (x,Q 2 ) [^-( n (Q 2 )-|)/n(l-Q 2 ) 

f» + (i-«i.(o))-^ 



(22) 



Now, if one takes W fixed with Q 2 ~ x — > 0, one can easily see that this equation simply 
reduces to: 

dF 2 (x,Q 2 ) 2 

^ ( 23 ) 

The calculations presented in the paper are based on the assumption of independent x 
and Q 2 , i.e. eqs.([20l, pT]) . In this context, the effective x-slope ^ e ff = d\nF 2 (x, Q 2 )/d\ia (1/x) 
is given by 



/x) ~ ™ >^ i- x vw )^-)\ (24) 



+ MO) - 1 + ^(Q 2 ) + B J B %_ X) ] ■ Fj 



NS- 



For Q 2 > Ql, the slopes have to be calculated from the evoluted structure function. 
Here, there are two fundamental procedures, the pure numerical and the mainly analytical 
calculations. The pure numerical procedure is very simple: 

dF ^Q2 2) « Q 2 ■ 2^ " mx, Q 2 + 5Q 2 ) - F 2 (x, Q 2 - SQ 2 )}, (25) 
dl d\lt/xP * (-^■ T ^^y 2 h-[Hx + 6x,Q 2 )-F 2 (x-6x,Q 2 )]. (26) 
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F 2 (x, Q 2 ) is the evoluted structure function whereas SQ 2 and Sx denote the corresponding 
increments which are fixed to be 10~ 3 • Q 2 or 10~ 3 • x in the calculations presented. For 
low Q 2 , we have checked this procedure comparing the values of eqs. (^,^) with those 
calculated using eqs. ( |20| , |2"1| ) . The agreement is very good which, in some cases, is demon- 
strated by identical numbers. This numerical procedure is the method used to determine 
the effective x-slope A e // = d\nF 2 (x, Q 2 )/d\n (1/x) of the evoluted structure function. 
In the case of d F 2 (x,Q 2 )/d In Q 2 there is, in addition, the way of mainly analytical cal- 
culations. If the parton distribution functions are known their derivatives concerning Q 2 
can be calculated from the DGLAP-equations. Instead of Q 2 the parameter 



S = ln{£}, T = \n(Q 2 /A 2 QCD ), T = \n(Q 2 /A 2 QCD ) 
is often used in perturbation theory. In terms of S 



dF 2 



dF 2 



d\nQ 2 \n(Q 2 /A 2 QCD ) dS 



(27) 



(2* 



and in the MS'-scheme 
ldF 2 (x, S) 



x 



dS 



dq(x, S) dq(x, S) 
dS dS 

ot s {Q 2 ) \ ^ fdq dq 



2ir 
1 da s (Q 
27 dS 



[C q ,2 *(q + q) + 2- C g , 2 * g] 



(29) 



The numerical integration procedure for solving the DGLAP-equations used in the 
work presented here gives the evoluted parton distributions and their derivatives on S as 
the output. In NLO, da s (Q 2 )/dS is simple to calculate, 



a s {T) 
2tt 
1 da s {T) 
dT 
1 da s 
2^~dS 



2tt 



2 / Aln(T)\ 
PoT\ ffi T J 
1 a s (T) _2A_ 
PIT* 



T 



T • 



2tt 

da* 



(ln(T) - 1) 



(30) 



2tt dT 



Thus, with the derivatives dq/dS, dq/dS and dg/dS one gets the Q 2 -derivative of F 2 . This 
method is called as "mainly analytical" (it includes a numerical integration procedure). 
Eq.flD 



11, 12 



i.e. 



Ql < Q 2 < Ql 



and in the 



is valid below the charm threshold 
high Q 2 -region where the charm can be considered as a "massless" dynamical quark [T^ 
As described above, the charm is treated via a photon-gluon fusion process in the range 
Q 2 < Q 2 < Q 2 = 50 GeV 2 [11], 0. From eq. ( |17D the charm slope contribution can be 
detrmined as 



1 dF%(x, Q 2 , m 2 c 



x 



c/lnQ 2 



2e 



a s (/i 2 
2tt 



1 dy 

ax y 



dC 



3,2 



x m: 



dlnQ 1 



Vf 2 )' g{y ^ 2) 



(31) 
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The total slope is the sum of eqs. ( [Zi^D and fl3"T|). 

We have calculated the Q 2 -slope of the evoluted F 2 in the perturbative region 
(Q 2 > Qq) using both, the numerical and the analytical methods. The values are in 
agreement although the difference increases somewhat in the region near to Qq. The 
values presented in the figures are from the numerical calculation. 
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Figure captions 



Figure 1. F 2 as a function of x computed in the CKMT model for twelve different 
values of Q 2 , and compared with the following experimental data (see ]l4j for the ex- 
perimental references): ZEUS SVX95 (black circles), HI SVX95 (white triangles), ZEUS 
BPC95 (white squares), E665 (white diamonds), and ZEUS 94 (white circles). The dotted 
line is the theoretical result obtained with the pure CKMT model, and the bold (solid) 
line is the result obtained with the NLO QCD-evoluted CKMT model when one takes 
Ql = 2.GeV 2 (Ql = A.GeV 2 ). 



Figure 2. F 2 as a function of x computed in the CKMT model for six (a) and five 
(b) different low values of Q 2 , and compared with the following experimental data (see 
HUH for the experimental references): ZEUS BPT97 (black circles), ZEUS BPC95 
(white circles), and E665 (white squares). The theoretical result has been obtained with 
the pure CKMT model. 



Figure 3. dF 2 / 'dlnQ 2 as a function of x computed by performing the NLO QCD 
perturbative evolution of the CKMT model (see appendices A and B for details on the 
calculation), and compared with the fit of the ZEUS F 2 data in bins of x to the form 
a + blnQ 2 (see reference [14| and references therein for more details on the data and the 



experimental fit). The dotted line is the theoretical result obtained with the pure CKMT 
model, and the bold (solid) line is the result obtained with the NLO QCD-evoluted CKMT 
model when one takes Q 2 = 2.GeV 2 (Q 2 = A.GeV 2 ). 

Figure 4. dlnF 2 /dln(l/x) as a function of Q 2 calculated by performing the NLO QCD 
evolution of the CKMT model, and compared to the fit F 2 = Ax~ Aeff of the ZEUS [|TJ[ 
and the E665 || data with x < 0.01. For details on the CKMT calculation, see Appendices 
A and B. The dotted line is the theoretical result obtained with the pure CKMT model, 
and the bold (solid) line is the result obtained with the NLO QCD-evoluted CKMT model 
when one takes Q 2 = 2.GeV 2 (Q 2 = A.GeV 2 ). 

Figure 5. dlnF 2 /dln(l/x) as a function of Q 2 calculated by performing the NLO 
QCD evolution of the CKMT model, and compared to the fit F 2 = Ax~ Aeff of the HI 
data (JT3J . For details on the CKMT calculation, see Appendices A and B. The dotted 
line is the theoretical result obtained with the pure CKMT model, and the bold (solid) 
line is the result obtained with the NLO QCD-evoluted CKMT model when one takes 
Q 2 = 2.GeV 2 {Ql = A.GeV 2 ). 

Figure 6. Gluon density distribution as a function of Q 2 calculated by performing 
the NLO QCD evolution of the CKMT model, and compared with the HI Collaboration 
data in reference [2D|. g(x, Q 2 ) is plotted and not xg(x, Q 2 ), in order to show more clearly 



the evolution with the scale. In the theoretical calculation, the bold (solid) line has been 
obtained by taking a value of Ql at the starting point of the QCD evolution, Ql = 2.GeV 2 
(Ql = A.GeV 2 ). 
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Figure 7. Gluon densities at /i 2 = 25.GeV 2 as a function of x calculated by performing 
the NLO QCD evolution of the CKMT model, and compared to those determined from 
HI DIS data (black dots) and from HI photoproduction data (stars). Experimental data 
on D* meson cross-section measurements are from references []T6L ^0[. In the theoretical 
calculation, the solid (dotted) line corresponds to a value of Qq at the starting point of 
the QCD evolution, Q 2 = 2.GeV 2 (Q 2 = A.GeV 2 ). 



Figure 8. Gluon density at fi 2 = 200. GeV 2 as a function of x calculated by performing 
the NLO QCD evolution of the CKMT model, to be compared with that obtained from 
the analysis of the HI di-jet data [16, 21 . In the theoretical calculation, the solid (dotted) 
line has been obtained by taking a value of at the starting point of the QCD evolution, 
Q 2 = 2.GeV 2 (Q 2 = A.GeV 2 ). 



Figure 9. Prediction of the behavior of xg(x, Q 2 ) as a function of x for several values 
of Q 2 measured both by HI and ZEUS collaborations. The experimental points are not 
shown since the analysis of the more recent data is not completed. The solid (dotted) lines 
have been obtained by taking a value of Ql at the starting point of the QCD evolution, 
Q 2 = 2.GeV 2 (Q 2 = A.GeV 2 ). 
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